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(yQ • We discuss some single particle properties of ^^O and ^'^Ca using correlated 



X 
J3 



basis function theory and Fermi hypernetted chain equations with central and 
tensor correlations. In particular, we concentrate on one body density matrix, 



> 

• momentum distribution, natural orbits and quasi hole states. The correlations 

o 



are variationally generated by a realistic hamiltonian containing the Argonne 
v'^ two-nucleon and Urbana IX three-nucleon interactions. The correlated 
momentum distributions show the well known enhancement at large momenta, 



o 

. 

2 ■ with a relative importance of the different correlations (Jastrow and tensor^ 

similar to that in nuclear matter. The natural orbits and their occupation 
numbers are obtained by diagonalization of the density matrix. The correlated 
first natural orbits occupations are depleted by more than 10%, whereas the 
first following ones are occupied by a few percent. The spectroscopic factors 
of the valence states are lowered by ~ 8 — 12% with respect to unity by central 
and tensor correlations, confirming that short range correlations alone are not 
able to explain the values extracted from (e, e'p) experiments. 
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I. INTRODUCTION 



The behavior of the single nucleoli in the medium bears clear signatures of the nucleon- 
nucleon (NN) correlations induced by the strong nuclear interaction. The most celebrated 
one is the presence of high momentum components in the momentum distribution (MD), 
which may be explained only in terms of short range correlations (SRC) and cannot be 
justified by any independent particle model (IPM) [|I|,0. In addition, both short- and long- 
range correlations (LRC) are now widely thought to be the origin of the large reduction of 
the spectroscopic strengths of the hole states 

In general, the investigation of SRC is one of the main topics in the physics of strongly 
interacting systems, from liquid Helium to nuclear matter. The translational invariance of 
homogeneous media has made possible a quantitative assesment of their importance in these 
contexts. However, the advances in realistic many-body theories have only recently allowed 
to address the same issue in inhomogeneous objects, as atomic nuclei, with a comparable ac- 
curacy. To this aim, correlated basis functions (CBF) theory has emerged as an effective tool 
to tackle the SRC problem in nuclear systems when realistic hamiltonians are employed. 
When used in conjuction with the Fermi hypernetted chain (FHNC) integral equations 
CBF has attained in ^^O and ^'^Ca the same accuracy as in the best variational studies of 
nuclear matter 

In order to dissentagle the SRC effects, one has to look at effects that go beyond any 
possible IPM description of the nucleus. SRC are mainly produced by the repulsive core of 
the NN interaction, whereas important intermediate and long range contributions (essen- 
tially of the tensor type) come from the exchange of one or more pions. This part of the 
nuclear interaction is the responsible of the nuclear binding . Direct evidence of the SRC 
are not easy to experimentally single out and only in these last few years, thanks to the 
advances in the experimental techniques, there have been consistent efforts aimed to their 
identification. For instance, (e, e'p) data in the quasi-elastic region need a large reduction of 
the IPM hole strength to be reproduced 0. Moreover, charge density distributions obtained 
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by elastic electron scattering experiments are smaller than those predicted by the IPM 



in the nuclear interior. These facts can be be explained by assuming occupation probabilities 
of the single particle levels different from the IPM ones [|Il| . 

The basic quantity to investigate in order to verify the hypothesis of partial occupation 
probability is the one body density matrix (OBDM), p(ri,ri'), defined as: 

p(ri,riO = {^o{A)\a^{rMri')\MA)), (1) 

where \l/o(^) is the ground state A-body wave function and a^(a)(ri) is the creation (an- 
nihilation) operator of a nucleon at the position ri. The Fourier transform of the OBDM 
gives the momentum distribution (MD), n{k), sometimes used in plane wave impulse ap- 
proximation to study inclusive and exclusive reactions. The natural orbits (NO) |12[, with 
their occupation numbers ( defined as the basis where the OBDM is diagonal. In 

the IPM, the nuclear ground state is described by a Slater determinant of fully occupied 
single particle (SP) wave functions below the Fermi surface, ap. So, the NO and the SP 
w.f. coincide, with na<ap=^ and ^a>ap=0. Deviations from this situation are a measure 
of the correlations, since they deplete the populated orbitals and allow higher NO to get 

Other important quantities are the quasihole (QH) wave functions, iph{'^)-, defined as 
the overlaps between \l/o and the hole states, obtained by removing a nucleon from 
the position r. From (e, e'p) experiments it is possible to obtain an accurate determination 
of the QH overlap functions |jl3[, whose normalizations give the spectroscopic factors, Sh- 



Typical values of Sh-, as extracted from the experiments, are Sh ~ 0.6 — 0.7 [0], and the 
deviations from unity (their IPM value) come in account of various effects, from center of 
mass to correlation corrections. 

In CBF theory the description of the correlations starts with an A-body correlated wave 
function 

vl>o(l,2...A) = 5(ni^.,)$o(l,2.../l), (2) 
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corresponding to a symmetrized product of non commuting two-body correlation operators, 
Fjj, acting on the mean field wave function, ^q{1,2...A), given by a Slater determinant of 
single particle wave functions, 0a (i). A realistic choice of Fij is: 

= E nr^J)0^, , (3) 
p=l,6 

where 

0^=''' = [1, a, ■ a,, (L ■ S),,] ® [1, r, ■ r,] (4) 

and Sij = (3 Vij ■ cTi iij ■ ctj — cri ■ <jj) is the tensor operator. 

The variational principle provides a natural recipe to determine the correlation functions, 
/^(r), and the single particle wave functions by minimizing the ground state energy. In our 
calculation we adopt a non relativistic nuclear hamiltonian of the form: 

^ = ^Ev.^ + E-. + E ^-^^ (5) 

i i<j i<j<k 

where we have used the v'g reduction of the Argonne fig ||15| potential and the Urbana IX 
|T6| three-nucleon interaction. 

II. ONE BODY DENSITY MATRIX, MOMENTUM DISTRIBUTION AND 

NATURAL ORBITS 

The one body density matrix (|T]) of doubly closed shell nuclei in (Is) coupling can be 
evaluated in FHNC [|T7|, starting from 

p(ri, nO = ^ I J d^A ^lH, 2, ..A)^,{1', 2, ..A) , (6) 

where Af = J d^ri...J (i^r^ |\l/oP and A is the mass number. Once the OBDM is known, its 
diagonal part provides the one body density, pi(ri), and the momentum distribution may 
be computed by 

n{k) = j I dhi I dh-v p(ri, nO e*k-(ri-r,,) _ 



The independent particle model expression for the OBDM is well known and reads as: 

P/PA/(ri,riO = E</'l(l)<^«(l')- (8) 

a 

In the case of Jastrow correlated wave functions (when only the first component of the 
correlation is retained), the density matrix may be expanded in powers of the dynamical 
correlations, h{r) = [f^ir)]"^ — 1 and uj{r) = f^{r) — 1, and of the statistical correlations. 
The expansion generates cluster terms classified according to the number of particles and to 
the number of the correlation lines. The FHNC equations allow for summing cluster terms 



at all orders. Details of the finite systems FHNC theory may be found in Ref. [^jIT^]. The 
FHNC equations for the more general /g correlation of were derived in Ref. |l7|,p!7[| for the 
one and two body densities and for the OBDM. 

The single particle wave functions needed to build the Slater determinant ^o{l,2...A) 
have been obtained by solving the single particle Schrodinger equation with a Woods-Saxon 
potential, 

^^^^^) = TT /° pw 1 ' (9) 

1 + exp [(r - Ro)/ao\ 

whose parameters have been fixed to reproduce at best the nuclei empirical densities. This 
choice, while spoiling the binding energies by only a few percent with respect to a full 
minimization, provides an accurate description of the nuclear densities §]. 

The correlated momentum distributions computed within the CBF scheme are shown 
in Fig.[l|. The MDs of ^^O (continuous line) and of ^°Ca (dashed line) are compared with 
that of nuclear matter with Jastrow (thin solid line) and full (dash-dotted line) correlations. 
The IPM momentum distributions are given by squares (^^O) and stars (^°Ca). We stress 
that the differences between the Jastrow and the /g correlations are similar in the infinite 
and finite systems and that the three cases show an analogous behavior in large momentum 
region, that is dominated by the short range structure (and the NN correlations) of the 
nuclear wave function. The effect is, to a large extent, independent on the nucleus. The 
tensor correlations enhance the tails of the MDs by a factor 2-3 with respect to the Jastrow 



case, slightly smaller than that found in Ref. |]T9[ (~4). The difference may be understood 
in terms of the stronger tensor force of the Argonne f 14 potential adopted in that reference, 
and of the presence of spin-orbit correlations, omitted in our calculation. 

Another piece of information that can be extracted from the OBDM are the natural 
orbits and their occupation numbers. The NO are obtained by diagonalizing the OBDM: 

Pi(ri,riO = E^«'/>r(ri)Vr(ri')- (10) 

a 

For spherical nuclei in {Is) single particle coupling, saturated in both spin and isospin, 
this expression can be recast as: 

2/ + 1 

pi(ri,ri/) = V— Pz(cos6'ii/)p;(ri,ri/) , (11) 

where Pi{x) are the Legendre polinomials, Oiy is the angle between ri and ry and z/=4 is 
the nucleon degeneracy. Exploiting again the spherical symmetry, one obtains 

pKri,riO = z^E^n.'/>^''(n)0r(riO. (12) 

n 

The first three NO of ^^O and '*°Ca are shown in Figs. ^ The NO occupation numbers 
for Jastrow and /e correlations are given in Table | The effect of the correlations on the 
shell model orbitals are mainly visible in the short range part of the Is state, making this 
NO more localized than its IPM counterpart. The shape of the other IPM states is barely 
influenced. The occupation of the NO corresponding to the shell model ones is depleted by 
as much as 22% (the 2s state in '^''Ca). In contrast, the mean field unoccupied states become 
sizeably populated. The two effects are largely due to the tensor correlations. 

The Green function (GF) approach of Ref. ||2^ found the ^^O n =1 NO more populated 



than the CBF ones for the occupied shell model states [l =0, 1), and, consequently, lower 
occupations for all the remaining orbitals (r;,^^=.921, ■^,^^=.941 and n^^=.017. The Ip {Id) 
occupation numbers are the average of the lpi/2 and lp3/2 (1(^3/2 and 1^5/2) orbitals of the 
reference). These discrepancies are probably due to the different potentials adopted, rather 
than to the methodologies. In fact, in the GF calculation the one-boson-exchange Bonn 
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B potential was used. The A8'+UIX model induces stronger correlations, so giving a 
larger depletion of the lowest NO. A similar effect was found in '^He atomic drops [^], where 
the strong repulsive interaction between the '^He atoms depletes the shell model occupations 
by 15-46 %. 



III. SINGLE PARTICLE OVERLAPS 

In a fixed center reference frame, the overlap between the A-body ground-state and the 
(A-l)-body hole state of the residual system, iph, is given by 

In CBF theory ^,^(1,2...^ - 1) is buih by substituting $0(1, 2. ..A) in (D with the Slater 
determinant obtained by removing a nucleon in the state h, ^h{l,2...A — 1). In doubly 
closed shell nuclei in the Is coupling scheme the radial dependence of the overlap function 
can be singled out, as 

i!hix) = i'hir)Yimif)xaT , (14) 

where Xcrr is the spin-isospin single particle wave function. In the IPM, the overlaps are 
simply the shell model functions and ipl^^^{r) = Rh=ni{f^), where Rh{fi) is the radial part 
of 

In order to develop a cluster expansion for il^h{f), its expression is rearranged as: 

Mr)=Mr)K^\ (15) 



where 



(16) 



and 



_ (vl/,(A-l)|vl/,(A-l)) 



Cluster expansions and FHNC-like equations are then used to compute Xh and Afh JT^ 



along the lines given in Ref. ||23[ to evaluate the overlap matrix elements in the CBF study 
of the nuclear matter spectral function. 

Finally, the spectroscopic factor, Sh, is given by the quasihole normalization, 

5, = 1 r'dri;l{r) . (18) 

The independent particle model gives Sj^^^ = 1. Center of mass (cm) corrections are a first 
source of correction. In fact, in the harmonic oscillator model they enhance Sh for the valence 
hole states (those having the largest oscillator quantum number, A^^,) by the [^/(yl — 1)]^" 
factor ||2^. So, the cm-corrected Ip-shell spectroscopic factor of ^^O is S^^^^^ = 16/15 ~ 
1.07, while the average between the 2s and Id states in ^^Ca is S^^i^^cm = (40/39)^ ~ 1.05. 



Similar results are found with the more realistic Woods-Saxon orbitals [|25[ 



The correlated spectroscopic factors (without cm corrections) in ^^O and ^''Ca are shown 
in Table |1| for /g and Jastrow correlations. Jastrow correlations marginally reduce Sh (at 
most 3%). Central spin-isospin correlations also provide a few percent depletion in the 
valence states, whereas the tensor correlations (/g) give most of the reduction. For instance, 
they lower Sip in ^^O to 0.90 and S2s and Sid in ^°Ca to 0.86 and 0.87, respectively. The 
Ip CBF ^^O result fully agrees with the variational Monte Carlo estimate of Ref. [^. The 
operatorial correlations largely influence the low lying states, whose spectroscopic factors 
are drastically reduced by both central and tensor components: Sis in is 0.70, Sip and 
Sis in ^°Ca are 0.58 and 0.55, respectively. 

The latest experimental values of Sp from the ^^0(e, e'p)^^N reaction p6[ are Sp^^^=0.61 
for the 1/2- ground state in ^^N and Sp3/2(6.32)=0.53 for the lowest 3/2^ state at 6.32 MeV. 
This state exhibits 87% of the total Sp^^^ strength, that is fragmented over three states at 
6.32, 9.93 and 10.70 MeV. The total Sp^^^ may be estimated to be ^^3/2=0.53/0.87=0. 61 



25| . In the Ca(e, e'p) K reaction ||T4| the transition to the 1^3/2 ground state gives 



Sd,^^ ~0.61±0.07. We recall that the CBF values are Sp =0.90 in ^''O and Sd =0.87 in ^^Ca. 
The squared overlap functions are shown in Fig.|^ for the Jastrow and /e correlations. 



and are compared with the IPM ones. Jastrow components have httle effect on the overlaps, 
while the spin-isospin correlations are the main responsible for the quenching of the IPM 
overlaps and of the spectroscopics factors. 

The IPM results obtained by a Woods-Saxon potential fitting the ^^0(e, e'p)^^N cross 
section to the 6.32 MeV state with S'p3^2(6.32)=0.53 are shown in the \ipip\'^-^^0 panel 
as stars. We have rescaled \ipip,FHNc\'^ by the factor 0.53/0.90 and the result is in nice 
agreement with the empirical estimate. 

IV. CONCLUSIONS 

In this contribution we have presented some results obtained within the correlated basis 
functions theory and concerning the properties of the single nucleon in the nuclear medium. 
More specifically, we have discussed the behavior of the one-body density matrix and some 
related quantities, as the momentum distributions and the natural orbits, in ^^O and *^°Ca 
using the Fermi hypernetted chain resummation technique. In addition, we have addressed 
the evaluation of the overlap functions and of the spectroscopic factors within the same 
approach. The relevance of the short range correlations, both of central and tensor type, 
has been stressed. The results presented here have been obtained using the Argonne v'^ 
two-nucleon potential plus the Urbana IX three-nucleon interaction, together with a set of 
single particle wave functions fixed to reproduce at best the empirical charge distributions 
of the two nuclei. 

The high momentum tail of the momentum distribution is largely dominated by corre- 
lations, and tensor components enhance the tail by a factor 3-4 with respect to the central 
ones. The tensor correlation is also important for the occupation of the natural orbits. In 
fact, the reduction of the occupation of the levels below the Fermi surface and, conversely, 
the enhancement above is amplified by the tensor terms. As far as the correlated overlap 
functions are concerned, they are close to the single particle wave functions if central corre- 
lations are used, whereas their shapes are strongly modified by the tensor correlations. The 
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overall spectroscopic factors are depleted by 10-15%, for the valence levels and 30-45% for 
the deeply lying ones. 

In spite of this reduction, the FHNC approach in ^^O does not reproduce the empirical 
^P3/2 spectroscopic factor extracted from (e, e'p) reactions. A similar situation was met in 
Ref. for nuclear matter, where the variational calculation of the one-hole strength, Z{e), 
around the Fermi level provided Zy{e ~ ep) ~0.88, mostly due to tensor correlations. Second 
order perturbative corrections in a correlated basis, including two-hole one-particle, {2h — 
Ip), correlated states, brought the strength to ZcBpi^ ~ ep) ~ 0.70. This decrease explains 
almost half of the discrepancy with the empirical ^°^Pb spectroscopic factor, Z(^^^Pb) ~ 
0.5 — 0.6. The missing strength can be attributed to the couphng of the single particle waves 
to the collective low-lying surface vibrations, not reproducible in infinite nuclear matter. It 
is expected that the inclusion of correlated 2h — lp corrections in the finite nuclei calculations 
can similarly take into account the coupling with surface vibrations. 
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TABLES 

TABLE I. Occupation numbers of the nl-th natural orbits for ^^O and ^°Ca in CBF, with the 
/e and Jastrow correlation models. 

nl nni{h-^^0) n„KJ;''0) n„K/6;'°Ca) n„KJ;'°Ca) 



Is 
2s 
3s 
4s 



0.858 
0.019 
0.010 
0.005 



0.960 
0.005 
0.002 
0.001 



0.864 
0.780 
0.052 
0.013 



0.952 
0.962 
0.002 
0.001 



Ip 
2p 
3p 



0.919 
0.021 
0.011 



0.980 
0.004 
0.003 



0.841 
0.024 
0.016 



0.949 
0.009 
0.006 



Id 
2d 

3d 



0.025 
0.011 
0.006 



0.006 
0.003 
0.001 



0.956 
0.030 
0.019 



0.983 
0.007 
0.006 



13 



TABLE II. CBF spectroscopic factors for ^^O and Ca, with Jastrow (J) and fe correlations, 
corr. Is Ip Id 2s 



J 0.98 0.98 

/e 0.70 0.90 

Ca J 0.98 0.99 0.97 0.98 

/e 0.55 0.58 0.87 0.86 
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FIGURES 
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FIG. 3. Squared quasihole wave functions. Solid lines: /e model; dot-dashed: Jastrow; 
dashed: IPM. The Ip panel of ^^O shows also the empirical overlap (stars) and the fe one, rescaled 
as explained in the text (lower dot-dashed line) 
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